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ABSTRACT 

Pulsars are rotating neutron stars, sweeping the emission regions from the 
magnetic poles across our line of sight. Isolated neutron stars lose angular mo- 
mentum through dipole radiation and (possibly) particle winds, hence they slow 
down extremely steadily, making them amongst the most reliable timing sources 
available. However, it is well known that younger pulsars can suffer glitches, 
when they suddenly deviate from their stable rotation period. On 2000 January 
16 (MJD 51559) the rate of pulsation from the Vela pulsar (B0833-45) showed 
such a fractional period change of 3.1 x 10~^, the largest recorded for this pulsar. 
The glitch was detected and reported by the Hobart radio telescope. The speedy 
announcement allowed the X-ray telescope, Chandra, and others, to make Target 
of Opportunity observations. The data placed an upper limit of 40 seconds for the 
transition time from the original to the new period. Four relaxation timescales 
are found, which are believed to be due to the variable coupling between the 
crust and the interior fluid. One is very short, about 60 seconds; the others have 
been previously reported and are 0.56, 3.33 and 19.1 days in length. 

Subject headings: stars: neutron - dense matter - pulsars: individual (PSR 
B0833-45) 



1. Introduction 

Observations of pulsar glitches, in addition to providing insights into the phenomenon 
itself, offer one of the few probes of neutron star structure, and thus the physics of ultra- 
dense matter. Vela is the brightest known radio pulsar, and as it is at a declination of -45", 



^School of Mathematics and Physics, University of Tasmania, GPO Box 252-21, Hobart, Tasmania 7001, 
Austraha 



- 2 - 



it is above the horizon at the Hobart Radio Observatory (Mount Pleasant) for more than 
18 hours a day. It undergoes large glitches in pulse rate every few years and so provides an 
excellent probe of the physics of glitches. 

Since 1981 the University of Tasmania has devoted a 14m diameter antenna at its Mt 
Pleasant Observatory to measurements of arrival times of pulses from the Vela pulsar. During 
this time we have observed 7 large glitches, or sudden decreases, in the period of the pulsar 
(McCuUoch et al. 1983, 1987, 1990; McCuUoch 1996). 

The telescope has a single pulse observing system whose speed and sensitivity have been 
enhanced in order to answer a number of questions; how quickly does the crust accelerate to 
the new period during a glitch, how soon does the recovery from the glitch start, and what 
is the form of this recovery? 

On 2000 January 16 (MJD 51559) the rate of pulsation jumped with a fractional period 
change of 3.1 x 10~^, the largest recorded for this pulsar. The glitch was automatically 
detected and we issued an lAU telegram (Dodson et al. 2000) within 12 hours, allowing 
the X-ray telescope, Chandra, to make Target of Opportunity (TOO) observations (Helfand 
et al. 2001). These observations have so far failed to find the signature of neutron star 
heating, which was the driver for the TOO, but have produced spectacular images of the 
X-ray pulsar wind nebula. 

2. Observations 

The reported observations were made simultaneously at 635 MHz, 990 MHz and 1390 MHz, 
to allow continuous measurement of the dispersion measure (DM). Dual channel receivers 
with an equivalent system temperature of 60 K are used at each frequency, except for 
1390 MHz, where only the right handed circular polarisation is observed. Receiver band- 
widths of 250 kHz at 635 MHz, 800 kHz at 990 MHz and 2 MHz at 1390 MHz hmit the pulse 
broadening from interstellar dispersion to less than 1% of the pulse period. The detected 
output from each receiver was folded for 2 minutes to give integrated pulse profiles from 
1,344 pulses. The profiles were subsequently combined to give profiles of total intensity. The 
signal-to-noise ratio in each total intensity profile is typically 30:1, allowing a mean pulse 
arrival time to be determined to an accuracy of 80/xs at 635 MHz, 60/xs at 990 MHz and 
180/is at 1390 MHz per integration. 

The recent improvements in the time resolution have been achieved by incoherently de- 
dispersing over 8 adjacent channels at 990 MHz, thereby increasing the signal to noise ratio 
and allowing observations of single pulses. The de-dispersed bandpass is sampled at 2 kHz 
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and recorded directly onto disk for later retrieval, while an 'on the fly' monitoring system 
folds on a ten second basis for which the RMS is 85/is. These arrival times are monitored 
and if a glitch is detected a warning is issued and the single pulse data are retained. 

3. Timing fits 

The full, folded, timing profiles are stored in EPN format (Lorimer et al. 1998), cross 
correlated with high signal to noise profile templates, then fitted to find the time of arrival 
(TOA) of the pulse. 

A Taylor expansion is parameterised in terms of the spin-down model phase (0) around 
a selected epoch {t relative to to) with frequency and frequency derivative z>, and has the 
form; 

(f>m{t) = 00 + Jy{t - to) + ^Z>(i - to)^ + ^i>{t - to)^ (1) 

The post glitch frequency change can be written as the sum of the pre-glitch frequency 
plus the permanent change at the time of the glitch, plus other decay terms: 

u{t) = i/o + Ai/p + Aupt + ^Ai>pt'^ + E„Ai/„e-*/^" (2) 

However, it should be noted that some investigators set the permanent change terms, 
Aup and Ai)p to zero. Wc have allowed Aup to take non-zero values, but keep Ai)p as zero. 
Both of these terms affect only the longest time-scales, where both timing noise and the 
occurrence of the subsequent glitches make definitive separation impossible. 

The data presented here were recorded between MJD 51505 and 51650. The arrival 
time data have been transformed to the Solar System baryccntre using standard techniques. 
The position and proper motion of the Vela pulsar was defined by data from the Radio VLBI 
position of Legge (2001). The recorded TOAs from all frequencies and both systems were 
fitted in the program TEMPO ^. The results of this fit are given in Table 1. 

Figure la shows the residuals from the pre-jump fit for data taken on 2000 January 16 
(MJD 51559) with the full polarisation system. Shortly after 07:34 UT, the residuals diverge 
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from the fit, indicating a sudden decrease in pulse period. Figure lb shows an hour of data 
starting at MJD 51559.3. Individual data points represent ten second averages constructed 
from the single pulse data. The period jump occurs on a very short timescale, without 
warning. The observations are consistent with an instantaneous change in period; modeling 
has shown that a spin-up timescale of forty seconds would produce a three sigma signal. 

The separation into four timc-scalcs is clear. The longer three decay terms are similar 
to those previously reported (Alpar et al. 1993; Flanagan 1990), and are in an approximately 
equal ratio of 5.9:5.7. These have been associated with the vortex creep models by Alpar 
ct al. (1993) and others. The fast decay timescale, not previously observed (or observable) 
is shown separated from the other effects in figure 2. We have subtracted the terms found 
by TEMPO in the 2 minute data from the single pulse data folded for 10 seconds. In this 
plot a gradual spin-up (as opposed to instantaneous) would be a negative excursion around 
the projected time of the ghtch, as we'd have overestimated the phase in the model. We see 
a positive excursion, followed by a rapid decay. A positive excursion could be produced by 
the pulsar slowing down just before the glitch or, more likely, if the estimate of the glitch 
epoch was too early because there was an extra component not resolvable in the 2 minute 
data. Once the original fit was removed this would give a linear rise with the gradient Az/ 
found in the 2 minutes data, followed by a decay. We have fitted this rise (AuAt) followed 
by a forth exponential decay term to give the true glitch epoch and the fastest decay term. 



4. The single pulse system 

Since the acceleration of the crust cannot be instantaneous, it should be possible to 
observe the spin-up of the rotation period of the pulsar. The parallel single pulse system 
designed to observe this has a three a detection limit of less than 40 seconds. This is 
in contrast to the observations made on the Crab (Wong et al. 2001), where the spin- up 
timescale has been observed to take about half a day. No variation could be seen either before 
or after the glitch in profile shape, magnitude, or polarisation in the radio observations on 
any time-scale. The raw time series using both long and short timescales has been closely 
examined to look for other periodic and quasi-periodic signals, but no extra terms were 
found. 
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5. Estimates for the force from the vortex lines 

We will use our data to explore one model of pulsar glitches in which the glitches are 
a consequence of pinning of superfluid vortices to the stellar crust. As the neutron star 
slows due the external drag on the crust, an angular momentum excess develops in the 
pinned superfluid. Ultimately the vortices unpin catastrophically, transferring the excess 
momentum to the crust, producing a glitch. 

Our high time resolution studies permit constraints on the manner in which the liquid 
interior imparts its angular momentum to the crust. This coupling is thought to be due to 
dissipative motion of vortices through the crustal lattice after unpinning. Various papers have 
tackled calculation of the drag on vortex hnes within the superfluid (e.g. Bildsten & Epstein 
(1989); Epstein & Baym (1992); Jones (1992)) produced by vortex-neutron scattering in the 
crust. We will ignore the details of the interaction here, as it is not a simple process. However, 
given that we have a lower limit for the spin-up time, we can obtain an observational limit 
for the force per length, /, that a vortex line exerts on the superfluid during this spin-up 
phase in Vela. Assuming that pinned vorticity in the crust is responsible for the glitch, it is 
straightforward to calculate the torque required to spin-up the pulsar, and then, by assuming 
canonical neutron star radii, crust depths and vortice lengths, give the force per unit length. 

During the spin-up, the torque on the crust, A/", is: 

N = IsAnjAt 

where Is is the moment of inertia of the portion of the superfluid that drives the glitch 
and Aflg is the change in angular velocity of that component. At is the spin-up time. Since 
angular momentum is conserved in the glitch, the torque is 

Af = I^ustAVt/At 

where Icrust is the moment of inertia of the solid crust plus any component that is 
coupled to it over a timescale in less than the spin-up time At, and AQ is the observed 
spin-up of the crust over At. 

In magnitude, the torque is 

where is the azimuthal component of the force that vortices exert on the lattice during 
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the glitch, and R is the average distance of vortices from the rotation axis, approximately 
the stellar radius. The azimuthal force is 

= fHN, 

where / is the force exerted per unit length per vortex, H is the average length of a 
vortex, and is the total number of vortices that are participating in the spin-up. 

The number of vortices is 

Ny = ny27rRAR 

where n„ is the number of vortices per unit area and AR is the thickness of the pinning 
region (approximately the thickness of the crust). The number density of vortices in a 
rotating superfluid is 

where fls is the superfluid rotational velocity, approximately equal to the rotational 
velocity of the crust and k is the quantum of circulation of a neutron fluid (0.0019 cm^/s). 
Putting these expressions together gives a constraint for /: 

AR, the thickness of the crust, is about 0.1 km. We take H, the length of a vortex, to 
be roughly AR. 

We defined Icmst as the moment of inertia of the solid plus any other component coupled 
to it in a timescale less than At. If the core fluid is tightly coupled to the solid crust over At, 
then Icrust is approximately the moment of inertia of the entire star (about 10'^^ g cm^). On 
the other hand, if the core is decoupled over At then Icmst would correspond to the moment 
of inertia of the solid crust only, which is about 10^"^ g cm^ for a moderate equation of state. 
So we can obtain a lower limit on /. 

/ > 10^^ IcrustA3 dyn/ cm 

Where I crust, 43 is the crust inertia in units of 10^^ g cm^. Vortices could pin to the lattice 
with a force ~ 10^ times this (Link & Cutler 2002), but our observations constrain the value 
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of / during the dissipative regime, whilst the pulsar is being spun up. Putting the same 
parameters into the Crab time-resolved glitch of MJD 50259.93 (Wong et al. 2001) we find 
that in that case / 4 x 10^ dyn/cm. This is further evidence, if more was needed, that the 
Crab pulsar is a very different object from the Vela pulsar. The major difference between the 
Crab and the Vela pulsars is the age, and therefore the temperature. Jones (1992) found a 
strong temperature dependence for some scattering processes, and these differences may be 
key to our understanding of the momentum transfer in glitches. Further spin-up observations 
for both the Crab and Vela are required. 



6. Future development 

Further improvements of the single pulse system arc being undertaken, with a new 
coherent dedispersion system, which will allow a detection limit of a few seconds, the order 
of the fastest coupling times in all theoretical models. This should also allow us to measure 
/, the force per unit length on the vortices, and confirm the fast decay term. Higher time 
resolution will allow further constraints on the coupling of the crust to the liquid interior, 
including the core. We will also co-ordinate with Chandra to reduce the TOO delay time in 
order to determine whether a very short lived thermal pulse might be emitted. 



The day to day maintenance of the 14m telescope, over 20 years has been a herculean 
task, to which too many people to name have contributed. Also we greatly appreciated 
the extremely helpful suggestions on the text from the referee, B. Link and also from S. 
Johnston. 
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Parameters for Epoch 51559 


u/Hz 


u/Hz s-^ 


i)/Hz 


11.194615396005 


-1.55615E-11 


1.028E-21 



Aup/Hz 


Aup/Hz 


3.45435(5)E-05 


-1.0482(2)E-13 


Tn 


Aun/lO-'^Hz 


1.2 ±0.2mins 


0.020(5) 


00.53(3) days 


0.31(2) 


03.29(3) days 


0.193(2) 


19.07(2) days 


0.2362(2) 


DM 


67.99 



Table 1: Parameters for the glitch epoch 51559.3190. Errors are the one sigma values. The data 
fit is from MJD 51505 through to 51650 (November 1999 to April 2000). 



MJD - 51559 



0.31 0.32 0.33 0.34 

MJD - 51559 



Fig. 1. — Arrival time residuals from the pre-glitch model for both the three frequency (la) 
and the single pulse (lb) systems. Figure la has folds of two minutes and figure lb has folds 
of ten seconds. 




Minutes 

Fig. 2. — The previously unobserved fast decay, with ten second folds. The other decay 
terms are removed revealing the later start, and decay of the fastest term. 



